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ABSTRACT 

We present a relatively simple operator formalism which reproduces the 
leading infrared logarithm of the one loop quantum gravitational correction 
to the fermion mode function on a locally de Sitter background. This rule 
may serve as the basis for an eventual stochastic formulation of quantum 
gravity during inflation. Such a formalism would not only effect a vast sim- 
plification in obtaining the leading powers of ln(a) at fixed loop orders, it 
would also permit us to sum the series of leading logarithms. A potentially 
important point is that our rule does not seem to be consistent with any 
simple infrared truncation of the fields. Our analysis also highlights the im- 
portance of spin as a gravitational interaction that persists even when kinetic 
energy has redshifted to zero. 
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1 Introduction 



Gravitons and massless, minimally coupled (MMC) scalars are unique in 
being massless without classical conformal invariance. The combination of 
these properties causes the accelerated expansion of spacetime during infla- 
tion to tear long wavelength virtual quanta out of the vacuum [HE]. As more 
and more gravitons and MMC scalars emerge from the vacuum, the metric 
and MMC scalar field strengths experience a slow growth. The effect can be 
felt by any quantum field theory which involves either the undifferentiated 
metric or an undifferentiated MMC scalar. 

An example is the one loop enhancement recently found [21 H] for the 
plane wave mode functions of massless, Dirac fermions which are coupled to 
quantum gravity on a locally de Sitter background, 

ds 2 = —dt 2 + a 2 (t)dx-dx where a(t) = e Ht . (1) 

(This background solves the classical Friedmann equation, (D — l)(d/a) 2 = 
A = (D — 1)H 2 , where A is the cosmological constant. The one-loop back- 
reaction on a(t) cannot affect the fermion mode function until two loop or- 
der.) At late times the full mode function H(ar;fc, s) behaves as if the tree 
order mode function E (x; k, s) was subject to a time-dependent field strength 
renormalization, 

{x;k,s) — > — . (2) 



This field strength renormalization takes the form, 

Z 2 (t) = 1 - ^-GH 2 ln(a) + 0(G 2 ) , (3) 

where G and H are the Newton and Hubble constants, respectively. 

The factor of ln(a) = Ht in expression ([3]) is known as an infrared log- 
arithm. Any quantum field theory which involves undifferentiated MMC 
scalars or metrics will show similar infrared logarithms in some of its Green's 
functions. They arise at one and two loop orders in the expectation value of 
the stress tensor and in the scalar self-mass-squared of a MMC scalar with a 
quartic self-coupling [5j. In scalar quantum electrodynamics they have been 
seen in the one loop vacuum polarization [6] and the two loop expectation 
values of scalar bilinears [7], the field strength bilinear and the stress tensor 
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[8]. In Yukawa theory they show up in the one loop fermion self-energy [H] 
and in the two loop coincident vertex function [TO]. In pure quantum grav- 
ity they occur in the one loop graviton self-energy [11] and in the two loop 
expectation value of the metric [12] . They even contaminate loop correc- 
tions to the power spectrum of cosmological perturbations [131 |H] and other 
fixed-momentum correlators [T5] . 

Infrared logarithms introduce a fascinating secular element into the usual, 
static results of quantum field theory. Their most intriguing property is 
their ability to compensate for powers of the loop counting parameter which 
suppress quantum loop effects. Indeed, the continued growth of ln(a) = Ht 
must eventually overwhelm the loop counting parameter, no matter how small 
it is. However, this does not necessarily mean that quantum loop effects 
become strong. The correct conclusion is rather that perturbation theory 
breaks down past a given point in time. One must employ a nonperturbative 
technique to follow what happens later. 

Certain models lend themselves to resummation schemes such as the 1/N 
expansion [IB] , but a more general technique is suggested by the form of the 
expansion for Z 2 (t) in 

Z 2 (t) = l + J2(GH 2 ) ^Q j0 [ln(a)] + Q, 1 [ln(a)] + . . . + c u _ x ln(a)j. (4) 

Here the constants are pure numbers which are assumed to be of or- 
der one. The term in (j3J) involving [GH 2 ln(a)] £ is the leading logarithm 
contribution at I loop order; the other terms are subdominant logarithms. 
Perturbation theory breaks down when ln(et) ~ 1/GH 2 , at which point the 
leading infrared logarithms at each loop order contribute numbers of order 
one. In contrast, the subleading logarithms are all suppressed by at least one 
factor of the small parameter GH 2 < 10~ 12 . So it makes sense to retain only 
the leading infrared logarithms, 

Z 2 (t)^l + J2ce,o[GH 2 \n(a)\ . (5) 

e=i 

This is known as the leading logarithm approximation. 

Starobinskii has developed a simple stochastic formalism [T7| which re- 
produces the leading infrared logarithms at each order [18] for any scalar 
potential model of the form, 

£ = -Idtfdvipg^yf^g - V(<p)V=jj ■ (6) 
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Probabilistic representations of inflationary cosmology have been much stud- 
ied in order to understand initial conditions [19J and global structure [20J. 
More recently they have been employed to study non-Gaussianity [21] . How- 
ever, we wish here to focus on Starobinskh's technique as a wonderfully simple 
way of recovering the most important secular effects of inflationary quantum 
field theory [22]. It is of particular importance for us that Starobinskh and 
Yokoyama have shown how to take the late time limit of the series of leading 
infrared logarithms whenever the potential V(ip) is bounded below [23]. This 
is the true analogue of what the renormalization group accomplishes in flat 
space quantum field theory and statistical mechanics. 

The solution of Starobinskh and Yokoyama is an amazing achievement, 
but it only gives us control over infrared logarithms which arise in scalar po- 
tential models ([6]). The most general theories which show infrared logarithms 
possess two complicating features: 

• Couplings to fields other than MMC scalars and gravitons; and 

• Interactions which involve differentiated MMC scalars and gravitons[j] 

An important step forward was the recent leading log solutions for MMC 
scalars which are either Yukawa-coupled to a massless, Dirac fermion [TO] or 
to electrodynamics [24J. Although the second model has derivative interac- 
tions, this feature was avoided (at leading logarithm order) by working in 
Lorentz gauge. We still do not understand how to treat derivative interac- 
tions. 

At the level of dimensionally regularized perturbation theory, the scalar 
leading logarithm solutions which have so far been obtained can be reduced 
to five simple steps [24J: 

1. Expand the full scalar operator <p(x) in powers of the free field (fo(x) 
which agrees with tp(x) and its first derivative at the beginning of in- 
flation as described in the recent paper by Musso [25] ; 

2. The expectation value of any desired operator can then be expressed 
as vertex integrations of retarded Green's functions times products of 
expectation values of pairs of free fields; 

1 Of course there would be no infrared logarithms if all the MMC scalars and gravitons 
were differentiated. However, infrared logarithms must arise, in the expectation values of 
some operators, from interactions which involve at least one undifferentiated MMC scalar 
or graviton. Examples include the h n dhdh interaction of pure quantum gravity I12j 
and scalar interactions of the form (p 2 dipdip [T3l I18j. 
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3. Make the following replacement for the expectation value of two free 
fields: 

O , H^T(D-l) 



<Po(x)<p {x') n) — ► - -g D 2 In min(a,a') ; (7) 

(47T) 2 L{^) 

4. Make the following replacement for the retarded Green's function: 

5. Evaluate the contributions from any other fields (for example, photons 
or fermions) exactly to the required order. 

Indeed, these rules even predict the occasional null results [26J that sometimes 
occur at low orders. 

It is straightforward to show that this old rule (fTHHJ) does not suffice to 
recover (jHJ). The purpose of this paper is to devise a simple rule which does 
work. We do not yet know if this rule applies either to other quantities 
or to higher loop orders. Nor do we possess a nonperturbative realization 
for this rule. Our rule nonetheless represents very significant progress in 
the struggle to solve inflationary quantum gravity in the leading logarithm 
approximation. Such a solution would make it simple to compute leading 
logarithm results at fixed order, and would also facilitate summation of the 
series of leading logarithms, thereby defining evolution past the breakdown 
of perturbation theory. 

In section 2 we explain how solving the Schwinger-Keldysh effective field 
equations is equivalent to computing the expectation value of a suitable 
canonical operator. Section 3 works out the operator and its expectation 
value to the order we require. At this stage the result is still exact and 
represents no simplification of the effective field equation technique. Our 
simplifying rule is presented in section 4. In section 5 we demonstrate that 
the rule indeed reproduces the leading infrared logarithm in the one loop 
correction to the fermion mode function. Our conclusions comprise section 
6. 



2 The Effective Mode Function 

We begin this section by describing the Schwinger-Keldysh formalism. This 
is a covariant extension of Feynman diagrams that produces true expectation 
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values instead of in-out matrix elements [271 EH1 ESI ED]- We then review the 
quantum-corrected Dirac equation whose solution (for spatial plane waves) 
gives the C-number effective fermion mode function Hj(x; k,s). The section 
closes by giving the connection between Hj(;r; k, s) and the fermion operator 

The in-out effective field equations give a fine representation of flat space 
scattering problems but they are not typically suitable for cosmological set- 
tings in which particle production precludes the in vacuum from evolving 
to the out vacuum. Persisting with thein-out formalism on de Sitter back- 
ground would result in processes being dominated by infrared divergences 
from the enormous spacetime volume of the infinite future [Tj [3T] • The bet- 
ter course in this case is to release the universe in a prepared state at finite 
time and let it evolve as it will. Problems of this sort are described by the 
Schwinger-Keldysh effective field equations [32] . 

Consider a scalar field (p whose Lagrangian (by which we mean the spatial 
integral of the Lagrangian density) at time t is L[(p(t)]. The fundamental re- 
lation between the canonical operator formalism and the Schwinger-Keldysh 
functional integral formalism is 



t(om)t{om)\$>) 



[d<p + ][d<p_]6 <pJfi)-<pJlti) ®*[(p_(t )]$Mt )} 



x0 2 [<p-]0i[(p+] exp 



dt {L[<p + (t)} ~ L[<P-(t)}} 



.(9) 



Here |$) is the Heisenberg state whose wave functional in terms of the p 
eigenkets at time to is <&[<£>(£o)]- The canonical expectation value on the 
left hand side consists of the product of an anti-time-ordered operator 2 [(p] 
times a time-ordered operator C?i [<£>]. The value of t\ > t is arbitrary as long 
as it is in the future of the latest operator occurring in either 0\ or (9 2 - 

The Feynman rules follow from relation (Q in close analogy to those for 
in-out matrix elements. Because the same field is represented by two different 
dummy functional variables, ip±(x) , the endpoints of lines carry a ± polarity. 
External lines associated with the anti-time-ordered operator 02 [y?] have the 
— polarity whereas those associated with the time-ordered operator Oi[cp] 
have the + polarity. Interaction vertices are either all + or all — . Vertices 
with + polarity are the same as in the usual Feynman rules whereas vertices 
with the — polarity have an additional minus sign. Propagators can be ++, 
— h H — or . 
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From this sketch we see that the N-point one-particle-irreducible (1PI) 
function T N (xi, . . . , x^) of the in-out formalism gives rise to 2 N different 
Schwinger-Keldysh 1PI functions T N (xi±, . . . , xn±)- Now recall that the in- 
out effective action is the generating functional of in-out 1PI functions, 

r M = E / ^kK^i) • • • / d 4 x N <p(x N ) x T N ( Xl , ...,x N ). (10) 

The analogous generating functional for Schwinger-Keldysh 1PI functions is, 

T[4> + ,(f)_] = / d^xx^xx) . . . J d 4 x N (p ± (x N )xT N \x 1± , . . . ,x N± ) . (11) 

The Schwinger-Keldysh effective field equations are obtained by varying this 
functional with respect to either <ft + or </>_, and then setting the two fields 
equal, 

StljJ =0 (12) 

It is worth being a little more explicit for the case in which the 0-point 
and 1-point functions vanish. If the classical action is S [</>], and the self- 
mass-squared is —iM± ± (x; x'), the Schwinger-Keldysh effective action has the 
following expansion, 

r[0 + , 0_] = sfa] - s[4l] -\Jd 4 xj dV 

<j) + (x)M* + (x; x')<j) + {x') + + (x)M + 2 _(x; x')(j)_(x') \ n( , 3 * , , 
+<P_{x)M 2 + {x; x')<f) + {x r ) + 0_(x)M_ 2 _(x; x')<f)_{x') J + W ] ' 1 ] 

The Schwinger-Keldysh effective field equations are, 

8S[<f>) 



- / dV 



M 2 (x; x') + M + 2 _(x; x') </>(x') + 0(</> 2 ) = . (14) 



The quantum-corrected Klein-Gordon equation results from linearizing (I14p . 
and its solution for a spatial plane wave is the scalar effective mode function. 
The peculiar combination of M^ + (x; x') +M 2 (x; x') in (fl4"|) has two important 
properties: 

• It is real, even though each self-mass-squared has a nonzero imaginary 
part; and 
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• It vanishes for any point x ,fl outside the past light-cone of x^. 

This paper concerns our solution of the quantum-corrected Dirac equation 
for the effective fermion mode function [4], 



x 



d A x' 



[x:x ) + 



(x;x'U S j(>) . 



(15) 



Hereby = 7^9 M and 7^ represents the usual, 4x4 gamma matrices. Note that 
Sj(x) is a 4-component (D-number field, even though the associated canonical 
operator ^i(x) is fermionic. There is no trace of the de Sitter geometry in 
the classical part of (|T5|) because we work in conformal coordinates, 



drj 2 + dx ■ dx = a 2 rj^ u dx^dx 1 ' where 



1 



-Ht 



Hr) 



(16) 



Massless fermions are conformally invariant in any dimension so we computed 
the fermion self-energy using dimensional regularization for the conformally 
rescaled field, 

%(x) = oT^ ijjiix) . (17) 



This removes any dependence upon the de Sitter scale factor from the tree 
order equation for ^(x), and hence for Sj(rc). 

Gravity is not conformally invariant, so one loop quantum gravitational 
corrections to the fermion self-energy involve the de Sitter scale factor. In 
computing these corrections we fixed the local Lorentz gauge so as to allow an 
algebraic expression for the vierbein in terms of the metric [34] • The general 
coordinate gauge was fixed to make the tensor structure of the graviton 
propagator decouple from its spacetime dependence [35]. After absorbing 
the divergences with three BPHZ (Bogoliubov-Parasiuk-Hepp-Zimmermann) 
counterterms we took the unregulated limit of D = 4 to obtain the following 
results [3J: 



in 2 H 2 ( ln(aa') 



K 



2 6 tt 2 [ H 2 aa> 
Tn(/z 2 Ax 2 H 



2 8 ttW 



Axl 



15 

— ln(aa/)<^— 71n(aa') 



>5 4 (x — x) 



+ 



k 2 H 2 f 15 2 -a 2 , 



, _ _ _x r ln(±F 2 A:r 2 j 

(-8^ 2 +4^V 2A 1 



Ax 2 + 



+ 7^V 2 



-ln(/i 2 Ax 



2 > 
++/ 



Ax 2 ++ 



Axl 



(18) 
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— K 



2%W Y 
+ (-8p 2 +4^V 



Ax 2 + _ 
-ln(±# 2 Ax 2 ; 



rln(/i 2 Ax 2 



Axl 



1 



Ax 2 _ 



Ax 2 + _ 



+ 0(« 4 ) . (19) 



Here k 2 = 167rG is the loop counting parameter of quantum gravity. The 
various differential and spinor-differential operators are, 



d 2 = rrd„d v , V 2 = dA , $ ee Yd, and 
The two conformal coordinate intervals are, 



1% 



Ax 2 + _(x; x) 



\x — x\\ 2 — I \T] — 7] 



i5 



\x- 



■x\\ 2 — (rj—rj' + i5^j . 



(20) 

(21) 
(22) 



Note that they agree for rj < rj', whereas they are complex conjugates of one 
another for rj > 7]'. 

Of course we can only solve for the one loop corrections to the field because 
we lack the higher loop contributions to the self-energy. Suppressing spinor 
indices and polarities, the general perturbative expansion takes the form, 



~(*) = 5> 2 ^) 



and 



£=l 



2/ 



(23) 



One substitutes these expansions into the effective Dirac equation (JT5 
and then segregates powers of k 2 , 







K Z I d A x' 



(24) 

Exl (x;ar , ) + [E 1 l (x; x')\ E (x') , (25) 



and so on. We considered the one loop correction Eu(x;k, s) to a spatial 
plane wave of helicity s, 

p-ikr) 



x\ k, s) 



Ui(k,s)e lk S where fc^f s) — k^Uj{k, s) . (26) 



In the limit of late times the source term on the right hand side of (1251) takes 
the form, 



k 2 H 2 
16tt 2 



K 2 i$Ei(x; k, s) 



17 

x —iHa'y E (x; k, s 



(27) 



Hence we conclude that the late time limit of the one loop correction to the 
effective mode function gives a time-dependent enhancement of the tree order 
field strength [3], 

{k 2 H 2 17 1 
l + j^-x—ln(a)\E (x;k,s). (28) 

We must now explain how the (D- number effective mode function 5(x; k, s) 
relates to the canonical fermion operator *f>(x). Consider the perturbative 
expansions of the Heisenberg operator equations for the graviton h^ v (x) and 
the (conformally rescaled) fermion ^i(x), 

h^(x) = h 0flI/ (x) + Kh lflv (x) + K 2 h 2 ^ u (x) + ... , (29) 
= ty i(x) + ^i l (x)+K 2 y 2i (x) + ... . (30) 

Long experience with such expansions permits us to anticipate how the first 
and second order corrections to \I/ depend upon the zeroth order fields, 

#i ~ h ^ , ^ 2 ~ h h ^ + # *o*o • (31) 

Because our state is released in free vacuum at t = (77 = — 1/H), it 
makes sense to express the zeroth order solutions in terms of the creation 
and annihilation operators of this free state, 

r d D ~ 1 k - 
h „u(x) = J ^o,! J2{ e ^(v, k, X)e lk ' s a(k, A) 

+e; u (r ] ;k,X)e' rk ^(k,X)}, (32) 



r cL k re ' -> -> 



v Vi (k,X)e- ik ^(k,s)} . (33) 



The graviton mode functions are proportional to Hankel functions whose 
precise specification we do not require. The Dirac wave functions Ui(k, s) and 
Vi(k, s) are precisely those of flat space by virtue of the conformal invariance 
of massless fermions. The canonically normalized creation and annihilation 
operators obey, 

a(k, A), a\k', A')] = 5 xxl {2i:) D - x 8 D -\k-k') , (34) 
{b(k, s), b\k\ s')) = 5 ss i27i) D - 1 5 D -\k-k') = [c(k, s), c\&, s')) . (35) 
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We can get the (D-number mode function Sj(x; k, s) from the zeroth order 
field ^oi(x) by anti-commuting with the fermion creation operator, 



ikr] 



E 0i (x;k,s) = {i&oi(x),tf(k,s)} = —=Ui(k,s)e 



ik-x 



(36) 



The higher order contributions to ^i(x) are no longer linear in the creation 
and annihilation operators, so anti-commuting the full solution ^%(x) with 
b'(k, s) produces an operator whose general form is, 



(37) 



The quantum-corrected fermion mode function we obtain by solving (JT5l) is 
the expectation value of this operator in the presence of the state which is 
free vacuum at t — 0, 



Ei(x;k,s) = (p\{^i(x),b f (k, s)} 



n 



(38) 



This is the promised relation between solving for the effective mode function 
and canonical operators [I]. 

Because we have a prediction (1271) for the late time limit of i $E(x; k, s) 
it makes sense to act the free kinetic operator on (J38l) . 



i$E(x; k, s) = {Q 



n 



(39) 



Of course this equation must hold order-by-order in the k expansions of 
E(x; k, s) and \fr(x). The order k° terms vanish identically. There is no order 
k, 1 correction to E(x; k, s), and the order k 1 correction to ty{x) vanishes when 
the expectation value is taken. The key relation for this paper comes from 
taking the late time limit at order k 2 , 



K 2 (n\{i@V 2 (x),tf(k,s)} 



n 



k 2 H 2 17 

—— x — iHa^Zoix; k, s) . (40) 
167r 8 



3 Perturbative Operator Solution 

The purpose of this section is to work out the canonical operator contri- 
butions to the left hand side of expression ( ]40l . We begin by giving the 
invariant action and fixing the gauge. This defines the fermion and gravtion 
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propagators which, in turn, give the retarded Green's functions. We then per- 
turbatively solve the Heisenberg operator equations to the required order in 
powers of the free fields f }32]|3"3l . Our result for K 2 i$ J &2 is reported in Table [TJ 
We also report the contribution of each term to K 2 {Q\{i jzSKE^^), b'(k, 
in Table [2j All the analysis of this section is done in D dimensions so that 
ultraviolet divergences are dimensionally regulated. 

The invariant Lagrangian density of Dirac + Einstein is, 

C = ^(^-(^-^^-^^^v^+^V^-^^^^v 7 ^ • (41) 

Here G is Newton's constant and H is the Hubble constant. The vierbein field 
is e^b and g^ u = e^ucV bc is the metric. The metric and vierbein-compatible 
connections are, 

i_lu = 2^ ^9o , fi,v 9va,ii 9hv,(t^ and A^ c d — c c(^ud,p ^^pu^pd 

(42) 

The Ricci scalar is, 

r^o^(t p -t p +t p f 7 -r p r CT 1 (A3) 

- LL — y ^ up,p 1 pp,i> ~ 1 pa 1 up x vu ppj ■ \* U J 

The gamma matrices 7^ have spinor indices i, j G {1, 2, 3, 4}, obey the usual 
anti-commutation relations and give the usual Lorentz generators, 

{ 7 6 ,7 c } = -2r/ fec I , J bc ^[ 7 b ,7 c ]. (44) 

It is useful to conformally rescale the vierbein by the de Sitter scale factor 
a(t), 

e pb = ae ph eP b = aT 1 & b . (45) 

Of course this implies a rescaled metric g^, 

2 " - a 2 {r]^ + nh^x)} where a = = e m . (46) 



9pv o g^ u 



The old connections can be expressed as follows in terms of the ones formed 
from the rescaled fields, 



IV = a-\6^a u +6^a,-g' } ' 7 a a g^)+V^ (47) 
A^d = -a~ 1 (e u c e lld -e u d e fic ^a^ + A^a . (48) 
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We define rescaled fermion fields as, 



^ = a 2 if) 



and 



2 if) . 



(49) 



We employ Lorentz symmetric gauge, = e^, which permits one to per- 
turbatively determine the vierbein in terms of the metric and their respective 
backgrounds [34] . 



J (3b 



(50) 



Here and throughout this paper graviton indices are raised and lowered with 
the Lorentz metric, e.g., h pu = r] pp ri ua hp a . The same convention applies as 
well to derivatives (d p = r]^ u d v ) and gamma matrices (7^ = 77^7^). The 
general coordinate freedom is fixed by adding the gauge fixing term, 

AC = -^a D - 2 V ^F^F u = ^(h ma ~h p ^+{D-2)Hah^ p 8f). (51) 

After some judicious partial integrations the gauge fixed Lagrangian density 
has the following expansion, 



K 



[1 h 2 --h pa h pa 

4 pa 



4 

1,, 1 . 1 1 



i ^p,. + -^ p ^ + i (^^p),. + i ^> w ,^7^ p ^ + 0(k 3 ) 

+-h^D fMU pa h pa + (Pure Gravity Interactions) . (52) 



The explicit form of the graviton kinetic operator D^J" 7 is not needed here; 
it can be found in ref. [3]. 

The ++ and H — fermion propagators are related to the conformal scalar 
propagator in the usual way, 



i[S] +± (x—x') = igUAf^x—x') =ify x 



rr£-r 



4vrf Axfr 2 



(53) 



The two conformal coordinate intervals Ax 2 +± were defined in f l2TH22l . 
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The graviton propagator takes the form of a sum of three scalar propa- 
gators times constant tensor factors 



pyAp a 



(x: x 



E l^iAUx-x') 



I=A,B,C 



(54) 



Because our gauge (J5l]) treats time and space differently it is useful to have 
expressions to the purely spatial parts of the Lorentz metric and the Kro- 
necker delta, 



V 



V = ^ + 5 /A and K = K- %S t 

With this convention, the three tensor factors in (15"4"1) are 

2 



rpA 
iiv L pa 

rpB 

V 1 pa 

rpC 

VV 1 pa 



2 V^pVaji 



D 3^ l|C " T 



2 



(55) 

(56) 
(57) 



(D-2)(D-3) 



P-3) ^ + V (^" 3 ) W + V • ( 58 ) 



0x0 



We follow the usual convention that parenthesized indices are symmetrized. 

The three scalar propagators in (1541) can be expressed in terms of the 
appropriate de Sitter invariant length function y+±(x; x'), 



y +± (x;x') = a(t)a(t')H 2 Axl ± {x;x') . 
The .B-type and C type propagators are hypergeometric functions, 

y+± 

4 

y+± 



iA* ± (x;x') 



H^TiD-mi) , D _ 21 D 
(4vt)t r(f) H 9' 



iA° ± (x;x') 



(4tt)t r(f) 
The A-type propagator has the intimidating expansion 
iA* ± {x- x') = iA ci 

h d-2 r ( jD _ 1 



2 

^rp-3)r(2) 2Fi(D _ 3 ^ ;l 



4 



(59) 

(60) 
(61) 



D 



(4tt)# r(f) \D-AT(D-l)^y +± 
lT(n+D-l),y +± 



2 — 7rcot(— d") + In (aa') 



fjD-2 oo 
+ — ~E 



(4tt) 



n=l 



[n r(n+f) V 4 



1 r(n+f + lW ± xn-f+2 



ra- 



il 



+ 2 T(n + 2) V 4 



(62) 
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We need retarded Green's functions in order to develop an expansion for 
the full fields in terms of the free fields of the initial time. There is a very 
simple relation between the retarded Green's function of any field and the 
corresponding ++ and H — propagators. If the field's kinetic operator is T> 
then the two propagators obey, 

ViA+ + (x;x') = id D (x — x) 



and T> iA + _(x; x' 
The associated retarded Green's function is, 







G(x; x') 



iA ++ (x; x') 



iA. 



From ( |63l) one easily sees that it obeys the required equation, 

VG(x;x') = -5 D {x-x') . 



(63) 
(64) 
(65) 



It also obeys the retarded condition of vanishing for rj < rf because the 
conformal coordinate intervals (I2T1) and (1221) are equal in that case. 
The Heisenberg operator equation for the fermion is, 



-K 



W| i ^--/ i >^- i (^M, CT - i ^>^]7^^vI/ + 0( K 3 ) . (66) 

We only require the analogous equation for the graviton to first order, and 
we only need the terms that involve fermions, 



D prj h 



K 

2 



OU 2 



+ (Pure Gravity Interactions) 
The next step is to expand Heisenberg operators in powers of k, 



h 



Hi/ 



fir 



(67) 



(68) 



I'Ohu + Khi pv + . . . . 

Of course the zeroth order equations (D^ u pa h 0pcT = and i$^o — 0) just give 
the zeroth order solutions h 0)Ml , and \& we already encountered in expressions 
and f)33p . respectively. The order k fermion equation implies, 



per 



pa 



*0 • 

(69) 
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Term 


Contribution to K 2 i^2{x) 


la 


-i/ t 2 /i ^(x) 7 ^^/d D a; / G c f(a:-a: / )C(^)7p^^o(a: / ) 


lb 




2a 




2b 


-i/ t 2 / io ^ Q (x) 7M J,^^/^x'af(^-^0C* /3 (^)7p^^o(x') 


3a 




3b 


± K 2 JW rG^](i;x')[^o(i')7p^$o(x')F x 7^^*o(^) 


4a 


-i/ ^ 2 9 a /^a; / rG'^](a;;a; / )^o(a: , )7p^^o(x , ) x ^J^o(x) 


4b 


\K 2 d a fd D x> rGP^( X] x')[%(x') % J aP M^ x 7/J ^o(z) 


5 


-InV(x) h p Q a {x)r llipll/ id^ {x) 


6 


-\K 2 h^ u {x) h^' a (x) r]^ a J ua ^ (x) 


7 


-±K 2 [/C(Z) C(^)]' Q ^p7,^a*0(^) 


8 





Table 1: Free Field Expansion of K 2 i^ 2 (x) 
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Hence the order k correction to the fermion operator is, 
%(x) = f ^x'G^x-x'^h^x')^ + CVhuJpJ^cK)- (70) 



In the same way we obtain the first order correction to the graviton, 
h1 v (x) - 



l - J d D xTG p °ix ] x')^ix') lp id'Mx') 



-d" 



^({x^lpJaa^x') \ + (Pure Gravity Terms) . (71) 



This brings us to the order k 2 correction to the fermion. We can of course 
drop any factors of i$*S?o — 0. With some further simplifications based on 
the first order equations we reach the form, 



3 



JV^+TtWT+k^lV^o • (72) 



Table[T]gives the free field expansion of i$^2, excepting only the contributions 
from the pure gravity corrections to hi^ v which vanish when the expectation 
value in PUj) is taken. 

Each contribution to Table [T] contains three free fields. It remains to 
evaluate the source term (HO]) . 



k 2 Q 



{i^ 2 ( x ),b + (k,s)}\n) + o(k a ) . 



(73) 



This is done by using the anti-commutator to absorb a \l/o and then exploiting 
the fundamental Schwinger-Keldysh relation (Q to express the expectation 
value of the two remaining free fields in terms of the propagator of appropriate 
polarity. To be definite, suppose the two remaining free fields are scalars 
(fo(x) and ip Q (x'). Here is where the factor ordering matters. From relation 
<Q we see that the H — propagator emerges from the order v^o^') X (p(x), 



n 



(p (x')(p (x) 



iA + _(x; x' 



(74) 



The order ipo(x) x ipo(x') gives the — h propagator, however, this is equivalent 
to the ++ propagator when account is taken of the factor of 6(rj—rj') in the 
retarded Green's function that is always present, 



G(x; x) x (il (p (x)(fo(x') 



n 



G(x;x') x iA_ + (x;x') , (75) 
G(x;x') x iA ¥+ {x;x') . (76) 
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Term 


Contribution to k 2 (VL [i^ 2 {x),b\k, s)} fl) 


la 


^ Jd D x'irA^] ++ (x; x>) ^d^G ci (x-x') lp d' a E (x') 


lb 


£ fd D x' d'Pi[^AP°} ++ (x; x') lil d v g)G c{ {x-x') lp J a pZ Q {x') 


2a 


^Jd D x' d a i[P v A ()CJ 


++ (x; x') 7 At J ra ^G ? c f(x-x')7 p ^S (x') 


2b 


f d D x > d a d'Hr^] ++ {x- x>) lp J va ^G ci {x-x') lp J a ^{x') 


3a 


if fd^x' [^G^]^; x') 7/ A^A; f _(x-x') 7p <9;:E (x') 


3b 


f idPxf df^G^^x; x') 7 M <%^Af_(x-x , ) 7p J (T/ 3So(x / ) 


4a 




4b 




5 


_3^i[^A^](x; x) 77 w 7^ ff S (x) 


6 


-f linv^ a to ir A" ff ] (x; x') ^ p7a J, CT S (x) 


7 




^A^](x;x)77 /ip 7,J (TQ So(x) 


8 


linv^ 9 to «rA^](x; x') r, palp J av Z {x) 



Table 2: Contribution to /t 2 (fi|{i^ 2 (a;), 6 f (fc, s)} 
free field expansion. 



f2 ) from each term in the 



17 



As an example we work out the (46) term. It is useful to begin by partially 
integrating the d' 13 without retaining the temporal surface term, 



(46) 



K 

T 



d D x' d a d' p 



(x] X') <&({x')j p J al 3^({x') X 7^ Jvc^ix) . 

(77) 

The term that contributes to the effective field equations is the expectation 
value of the anti-commutator of (46) with tf(k, s), 

.2 



{(46), 6^,*)} 



K 

T 



d D x'd a d' p 



(x; x 1 ) 



Hx'hpJ*p{*&'), b\k, s)} x lfM J ua ^x)\n) , (78) 



— j d D x'd a d'P[ lu G pa ](x;x') 



^ix')l p Jop 2d>') x ^J va ty({x) 



(79) 



At this stage the spinor indices become confusing so we write them out 
explicitly, and also remove all C-numbers from the expectation value, 

.2 



n 



Ab) h b\k, s )}\n) = ^ j d D x'd a d ,f3 



00 - 00 J 



K 

T 



d D x'd a & fi 



(x;x') x (y^Jva),, x {-f p J a ^ k E e(x') 



(1U fip<J 



%k(x')%j(x) 



n 



(80) 



(81) 



The expectation value on the final line of (IHTj) is minus the H — fermion 
propagator, 

(x;x') = —i$j k iA'f{x — x') . 



n 



%kix')%j{x 



n) = -i 



(82) 



The minus sign derives from the fact that the preferred order for the fermion 
propagator is Substituting (!82|) into (IHTI) gives an expression we can 

write without resort to explicit spinor indices, 



{(46),6t(M}|0) 



I d L> , ■'<)"■()" 



x;x')^ ti J 1/a ^iAllx-x')-f p J a pE (x') . (83) 



Table [2] gives our results for each entry in Table [U 
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4 Our Rule 



The first eight entries of Table [2] provide a somewhat cumbersome re-expres- 
sion of the nonlocal contributions to the order k 2 source term of expression 
( I25p . The original source has the generic form of a difference of ++ and H — 
terms, with each polarity being a product of contributions from the graviton 
and contributions from the fermion. Table [2] effects the following re-grouping, 



From expression we see that the difference of ++ and H — propagators 
for any field gives —% times the retrarded Green's function of that same 
field. The first eight entries come in pairs of this form: (la)-(3a), (16)-(36), 
(2a)-(4a) and (26)-(46). This is an illuminating insight but it represents no 
simplification of the original calculation. 

We cannot simplify the propagators and retarded Green's functions asso- 
ciated with the fermion. In contradistiction to the graviton, the fermion is 
a "passive" field which cannot produce infrared logarithms [TOj [21] . Passive 
fields contribute factors of order one that derive from both the infrared and 
the ultraviolet. To correctly recover these factors the passive field must be 
treated exactly. 

Our simplification concerns the propagators and retarded Green's func- 
tions of the graviton. The A-type graviton polarizations are the "active" 
fields which cause infrared logarithms, whereas the -B-type and C-type po- 
larizations are passive. Because this particular calculation involves only one 
graviton propagator or Green's function there is no chance of getting an 
infrared logarithm unless the v4-type part of the graviton propagator is in- 
volved. 

Even within the A-type polarization, only the following tiny portion of 
the infinite series expansion (|62|) of «A^ ± plays any role in generating infrared 
logarithms, 



/ 



d D x' {{ ++ ) h x - (+-), x ( + _)^} Hq^') = J d D x' (++), 

x{( ++ )^-( + -)^}H (x')+ Jd D x'{( ++ ) h -( + -) h }( + -UZ (x>). (84) 



i5A+±(x; x') 



h d ~ 2 r(p-i) 
(47r)f r(f) 




(85) 
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Our rule is accordingly to make the following simplifications on the graviton 
propagators and Green's functions, 



pumper 



^Tf xi6A*(x;x') 



x SGaIx; x') 



where the A-type tensor factor is (|56|) and we define SGa(x; x') to be, 



5Ga[x] x') 



idA^x; x') - i8A+_(x; x') 



1H 2 r(f +i) 



57T 2 



¥ D-A 



(aa!) 2 -% 



(aa'f 



AxPr 4 Ax?- 4 



(86) 
(87) 

(88) 
(89) 



In the next section we demonstrate that applying our replacements ( IHoTIHTl) 
to the various terms in Table [2] reproduces the source term (j4"0)) whose inte- 
gration gives the infrared logarithm ([3]). 

We close this section by commenting on the relation between our rule 
and the replacements ([THH1) that have been shown to reproduce the leading 
infrared logarithms to all orders in scalar models without derivative couplings 
[TO"! |23| [21]. The rules are certainly not identical but they do seem to agree, 
at leading logarithm order, for D = 4 and for certain treatments of the spatial 
coordinate separation. To see this, first take the D = 4 limits of (1851) and 



lim idAiJx; x') 



lim SGa(x; x') 

D^4 



H 2 



8tt 2 



In 



l^Ax 2 +± 



H 2 

- — 9[r] — r]' — \\x—x'\ 



If we set x! = x in fl90l) the result is 



(90) 
(91) 



IP 

8^2 



In 



4W 



1\2 

a 



+ 



£1 

471-2 



hi 



mm a, a 



+ 0(1) 



(92) 



At leading logarithm order this indeed agrees with the D = 4 limit of our 
previous rule (171). Similarly, the spatial integral of ( 19~TI) is, 



/ 



H 2 

d 3 x - — 9(7]— n'- 



\X — X 



1 



3#W a> 3Ha' 3 



1 + 0$ 

\fl+ 



(93) 
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At leading logarithm order this agrees with the D = 4 limit of the spatial 
integral of 

These correspondences seem to mean that our new replacements f l86ll87j) 
would reproduce the leading infrared logarithms of the simple scalar models 
previously studied. However, it is straightforward to check that the old 
replacements ([THED do not reproduce the result (T4U1) we get from quantum 
gravity, whereas our new replacements (18611871) do. It therefore seems that 
our new rule represents a successful generalization of the old rule to the 
more singular environment that arises when derivative couplings are present. 
What is not yet clear is whether or not the rule can be simplified. 



5 Analysis 

The purpose of this section is to show that applying our rule (I86H87P to Ta- 
ble |2] reproduces the result fj40l) of our explicit computation. We begin by 
observing that any terms involving derivatives of So cannot contribute at 
leading order. That reduces the problem to considering the nonlocal contri- 
butions 16, 26, 36 and 46, and the local contributions 6, 7 and 8. The local 
contributions were evaluated in an earlier effort to understand our result 
(|40|) on a qualitative level by making the Hartree approximation [4], so we 
concentrate on the nonlocal contributions. We first introduce a systematic 
classification for the myriads of distinct terms they give when the A-type 
tensor factor and the factors of 7 M J VQ . and ^pJ a fi ar e broken up. Then we 
explicitly evaluate four of the contributions from (26) as an example. Final 
results for all nonlocal and local contributions are reported in tables. 

It is important to understand that we only seek the leading late time 
behaviors of the various source terms in Table [2j By considering the form of 
quantum gravity interactions we see that the one loop mode function can be 
enhanced by at most a single infrared logarithm [24J , 

K 2 (tt\{<I?2(x),h\k,s)}\n) ~ k 2 H 2 x ln(a) x E (x) . (94) 

The source terms of Table [2] should be % ft times this, which gives the loop 
counting parameter k 2 H 2 times iaH'y 'E,o(x) . 

Now consider how the derivatives of Table [2] act. Any which act on the 
tree order mode function Ho(x') will bring down factors of the wave number 
k. This factor of k must persist, even after the integration over x' 11 , because 
the integral remains finite for k = 0. Further, this factor of k will always be 
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accompanied by a factor of 1/a to make the wave number physical. It follows 
that the fastest growth possible for any d'S (x f ) term is ln(a)/c7°S (a;). We 
can therefore forget about nonlocal contributions from (la), (2a), (3a) or 
(4a), and also the local contribution from (5). For the same reason we can 
make the following simplification in the nonlocal contributions from (16), 
(26), (36) and (46), 



E (x'; k, s) = E (x; k, s) x e 



E (x; k, s) x 1 . (95) 



We turn now to the problem of classifying the many distinct contributions 
that derive from (16), (26), (36) and (46). These four terms all involve a single 
factor of the A-type tensor and either one or two factors of the Lorentz 
generators. The A-type tensor indices are purely spatial, for example, 



X (liJja) X (jfcJ^j) 



(96) 



Our classification system is based upon decomposing the 7 • J factors as 
follows, 



1 1 1 



(97) 



The totally anti-symmetrized term drops out because the A-type tensor fac- 
tor is symmetric in i and j. We label the other two terms by Roman numerals 
"I" and "II" as follows, 



I 

II 



2 $ialj 



and 
and 



2^/37£ , 



(98) 
(99) 



The indices a — which appears only in (26) and (46) — and (3, contract 
into derivative operators <9 a and d' 13 that act upon the graviton propagator 
or retarded Green's function. The indices on I-type terms must be spatial 
but those on II-type terms can be either spatial — denoted by just II — or 
temporal — denoted by II'. The type-II term always produces a contraction 
of the A-type tensor factor, for example, 



l]rpK 



x 5 k £ = 



D-3 



5 13 . 



(100) 
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However, the type-I term can receive distinct contributions from each of the 
three terms in [ y T^]. Where the results are distinct we label these "A", "B" 
and "C" as follows, 



A 



B 



C 



D-3 



5 ij 5 



(101) 



The various classifications are arranged in a prescribed order, and are 
separated by periods. First comes the term designation — lb, 2b, 3b or 4b. 
Next comes the leftmost of the 7 - J factor designations — I, II or II'. If there 
is a second 7 • J factor, its designator comes next. The final designator is the 
A, B, or C from the tensor factor, with no designator denoting the presence 
of all three terms. As an example, consider the full (lb) term, 



K 



— J d D x> d ,(i i5A? + (x- x') ^T^d^G^x-x'^J^x 1 ) 
The lb. II' contribution is, 



(102) 



lb.ir = 



K 



J d D x' x (&°) x i5A? + (x; x') x (-73-^ 

i 
2 



X' 



'^idjQG^x-x') x (--70) x So(x') , (103) 
2{ J^Jd D x'd' z5A^x ] x')WG^-^h ^o(x') . (104) 



IK 



In contrast, the contribution from Ib.I.B is, 

.2 



Ib.I.B = J d D x' x (& k ) x i5A? + (x; x') x (5 u 5 jk ) 

x-fidjfiGtfix-x') x Q 7< ) x E (x') , (105) 

= -^Jd D x'd k i5A^x;x'h%$G ci (x-x')YZ (x') . (106) 

Note the minus sign from converting d' k to —d k . 

To describe the evaluation technique we have chosen four of the contri- 
butions from (26): 2b.II'.I, 2b.II.I, 2bT.IT and 2b.I.II. Each of these involves 
a single derivative of the conformal Green's function, 



QG ci {x-x')- 12 



_ D 

2n~ 



+ 



Ax? + Ax?_ 



^A Xfl . 



(107) 
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Term 


Contribution from (n {(lb), tf(k, s)} ft) 


Coef. 


lb.II' 


J ^fd D x'd' iSA^ + (x; x')@$ G c{ (x -x')^Eo(x') 


+\ 


lb.II 


-^^I^x'dkid^x-x'^G^x-x')^^') 


1 

4 


lb.I.A 


— lj- / d D x'd k i5A+ + (x] x')^ k d^G c{ {x-x')^'^Q(x') 


1 

16 


lb.I.B 


-fjd D x'd k t5A^ + (x;x') 1 %^G c{ (x-x') 1 ^o(x') 


3 
16 


lb.I.C 


I ^Jd D x'd k i5A^ + (x-x')^G c i(x-x') 1 k E (x') 


+1 


Total 




1 

8 



Table 3: The full contribution from each (16) term consists of its numerical 
coefficient times ! j^r 



They also involve two derivatives of the A-type propagator, 



d k d' Q i5A^ ± (x;x') = 
d k d i5A+ ± (x; x) 



in 2 ( aa') 2 2 



f (D-2)Ar] (D—A)Ha! 



2Ax°- 2 



87r 2 (aa') 2 

dtdtmUvJ) = — U:L - 



(£>-2)Ary (D-A)Ha \ 
2A.r" I 



Ax? ± 



5 



kt 



Ax?~ 2 



+ 



[D-2)Ax k Ax k 
Mi 



(108) 
(109) 
(110) 



87T 2 (aa') 2 2 

In these and all subsequent expressions we define the coordinate differences, 

Ax^ = x^-x'^ and Ar] = 77 - rf . (Ill) 

Each of the four terms we are considering takes the form of a common 
integral operator acting upon a different integrand. The integral operator is, 



K 2 # 2 r(f)r(f +1) 



d D x'{aa'f-^ 



6Att d D-3 
The four different integrands are, 

2b.II'.I 7 7 V2o(a;') x A Xll Ax k < 

24 



1 1 

+ 



Ax° + Ax®_ 



x 



(D-2)Arj (D—A)Ha\ 
+ 2Ax°~ 2 j 



Ax» + 



(112) 



, (H3) 



Term 


Contribution from (to {(3&),&t(jfe, s)} fl) 


Coef. 


3b.II' 


2 ™\ ) Sd D x'& G dG A {x-x')$g)iA% (x x')^°E (x') 


+ 1 


3b.II 


-2(fe) IdPJdkdGAx; x'^iAf^x-x'^Eoix') 


1 

4 


3b.I.A 


-if Jd D x'd k dG A (x; x'^d^iAfX^-x'^E^x') 


1 

16 


3b.I.B 


-f fd D x'd k SG A (x; x')Yd k $iAi_(x-x')YZo(x') 


3 
16 


3b.I.C 


Id D x'd k SG A (x; x'W^tAf^x-x'^Eoix') 


+1 


Total 


^/d D x / ( 9 / ^ix;xO[ ^J T^]7 ^ 9#A c + f _(x-a; / )7fc^So(a: , ) 


1 

8 



Table 4: The full contribution from each (36) term consists of its numerical 
coefficient times lJ j^i 



2b .ii.i =, /yys„M x AxJ-^j + ( ° • (H4) 



2b ,. ir ^ 7V/ - M x ^{.(^ _(|«}, (115 ) 

2U.II => 7 V T <2„M x A % {-^ + <g "^ A / A3: ' l ■ (H6) 



If we ignore the difference between x' M and x M in the wavefunction E (x') and 
perform the angular averages, the various integrands take the form, 



2i,n.i = i- (u8) 



7 °S (x) x < 


f (D-2)A77||Af || 2 




7°So(*) x < 


f(D-l)A?7 (£> 


[ Ax°~ 2 


7°H (^) x < 


' (D-2)Ari\\Ax\\ 2 


[ Ax° + 


7 °S (a;) x < 


'(D-l)Arj (D 


[ Ax*' 2 



2b.II.I =^ 7 ^o(^) x <; ^ 2Ax D ~ 2 I ' ^ m > 



2b ,,r ^ )x r ^-" + ( °-^gf £ " 2 }. (us) 

2b.I, II /=„(*) x < _ g -gll^f } . (120) 

Each of these four terms can be written as a common factor times a sum 
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of integrals. The common factor is, 

'•• i // 2 r(f)r(f+i) _ 



64ttt + 2 D - 3 
The four fundamental integrals are, 



h 
h 
h 
h 



(D-2) J d D x' (aa') 2 ~% 
^(D-A) J d D x' (aa') 2 -T 
(D-l) f d D x' (aa') 2 ~^ 
l -(D-A) J d D x' (aa') 2 ~ 



7°So(^) • 



1 1 

+ 



Ax° + Ax»_ 

1 1 

+ 



Ax° + Ax®_ 

1 1 

+ 



Ax° + Ax®_ 



Ax° + Ax°_ 





(121) 


Ar]\\Ax\\ 2 

Ax° + ' 


(122) 


" Ha\\Ax\\ 2 


(123) 


Ax ++ 2 ' 


Ar] 
Ax°- 2 ' 


(124) 


" Ha'\\Ax\\ 2 


(125) 


Ax°- 2 



And the four terms under consideration are, 

.2 tj2 rfD\ r (D 



2b.II'.I 
2b.II.I 
2b.I.lT 
2b.I.II 



(126) 
(127) 



^ 2 £(|)r(f+i) f , 

* 647T* D-3 7H,(l)X ( Jl " /2 t' 

♦ ^S^**) *{*-',}. 

The procedure for evaluating Ji_ 4 is, 

1. Perform the angular integrations; 

2. Note that (for 5^0) the radial integrand vanishes for r > At]; 

3. Make the change of variable r = Ar^^/a 7 , which reduces the radial inte- 
grals to beta functions; and 

4. Make the change of variable rj' = —l/(Hat). 
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As an example, consider I\. The first step brings it to, 



h = (D-2) x x / d V '(aa') 2 -- / 



drr 



D-2 



1 1 ] Arir 2 

+ r o 71 . - ' • (130) 



i [r 2 — (Ar] — z5) 2 ] 2 [r 2 — (Ar/+i5) 2 ] 2 J [r 2 — (Ar] — z<5) 2 ] 2 

Step 2 is accomplished by noting that the =pz<5 factors serve to fix the phase 
of the complex numbers that must be raised to the D/2 power on the final 
line of f)130p . For r > Ar] that phase is zero, whereas it is ±7r for < r < A77, 



r 2 - (At; =f i$) = e x [Ar] 2 - r 2 j for < r < Ar] . (131) 
Hence the curly bracketed term of (11301) becomes, 

1 1 2isin(^) 

I ■ w = ^-4r for < r < Ar] . 



[r 2 — (Ar] — i5) 2 ] 2 [r 2 — (Arj+iS) 2 ] 2 [Arj 2 — r 2 ] 2 
It follows that steps 2-4 give, 

27f 2 /"J? _D 

Ji = (D^x—^x^dv'iaa'f—Av 



(132) 



/■Ar? n 2isin(^)e-^- , , 



< D - 2 > x THEE x />'K) 2 ~ f AT D+2 

dnf-^)e- i7r T fax x 2 * 1 (1 - x)~ D , (134) 
v 2 7 jo 



x« sin 



(D-2) x x H^ajldttf-^l-t)-™ 



r(¥) 

/Dttx . d r(^±i)r(-D+i) 

xirin(— )e— t x 2 ' I > (135) 



2 7 r( 

After step four the other three integrals are 



D-l 



1(73-4) x x H D -*a£dt#-*(l-t)-™ 



re_ 

27 



xirin(— )e-^T-Dx l2 r( U +5) ^, (136) 



H^>-^ >< ^IS"^ (137) 



2 / r(- 



x iS in(-)e-^-D x g_L^ (138) 



2 / r(- 

Note that (for D = 4) the £ integrands of Ji, J3 and J4 are finite at £ = 0. 
This means we make only an error of order 1/a by extending the range of £ 
down to £ = 0, at which point we get another beta function, 

jfa^fl-*)-™ = r(f -% +3) + °0- (139 » 

f^-')-- r( V-iT 4) + <)- (140) 

This allows us to evaluate Ji, J 3 and J 4 . Setting D = 4 — e and taking e to 
zero gives the following results for these three integrals, 

Ji — ► 2 x 47r x aH x — - x — i— e x — = i^ 2 x mi? , (141) 

2 2 16e 8 > k J 

X 7T X 3 

J 3 — > 3 x 47r x aH x — x i-e x = -7r 2 x iaH , (142) 

2 2 8e 8 v ; 

X X vi 3 3 

L — > — e x 47r x aH x - x z— e x — = — tt 2 x iaH . (143) 
2 e 2 8e 8 K J 

This procedure is no£ valid for the £ integral of I2 because the integrand 
diverges at £ = 0. The right way to evaluate the £ integral in (11361) is to first 
add and subtract the £ integral from (j!38j) . 

/ 1 ci££f- 3 (l-£)-^+ 3 

a 

= j'dtt%- 2 (l-t)- D+3 + jldt{—}t%- 2 {l-t)- D +*. (144) 

a a 

28 



Term 


Contribution from (to {(2b),bUk, s)\ Si) 


Coef. 


2b. II'. II' 


in 2 (D-l) ffjD /a ai-s:AA ( r ')^0an ( a 0^ ( a 


+0 


2b.ir.II 


fd D x'd n dki5A^Jx- x'WdGrdx x'^Znix') 


3 
4 


2b TT TT' 

_ YJ . A.A.. LA. 


« 2 (D-1) [jDlXfl/jSjAA ( r ')^ykgn ( i\ ow n ( r i) 
4(_D— 3) J i^fe^o'" ++v ' / ' V^^cf^ x / / "(A x J 


+ 8 


2b.II.II 


^^/A' ( 9 A; ^A4 4 + ( a ;; a :')7^G' cf ( a ;- a ;')yS ( a ;') 


3 
8 


2b.H'.I 


-^ ) Sd D x'd Q d k i5A^ + {x-x') 1 ^G c{ {x-x') 1 k ^{x') 


+t 


2b.II.I 


- I ^Id D x > d k d e iSA^ + (x;x') 1 k ^G cl (x-x')^o(x') 




2b.I.II' 


I ^jd D x'd k d' tSA^ + (x ] x') 1 k ^G ci (x-x') 1 °E (x') 


l 

8 


2b.I.II 


-^^ ) Sd D x'd k d,i5A^ + {x-x') 1 k ^G ci {x-x') 1 ^{x') 


+ 1 


2b.I.I.A 


-^jd D x'd k d k i5A^x-x') 1 ^G ci {x-x') 1 ^ {x') 


~32 


2b.I.I.BC 


-'-^E^I^x'dAidA^x-x'^^G^-x')^^ 


1 

32 


Total 


-i£fd D x'd a d'Pi8A£ + (x- x') [^T k %J jc ^G ci (x-x') lk J ep E (x') 


5 

16 



Table 5: The full contribution from each (26) term consists of its numerical 
coefficient times %l fJ^ Ha^E^x). 
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Term 


Contribution from {Q [(46), tf(k, s)\ fi) 


Coef. 


4b.IF.ir 




+0 


4b.ir.II 


^^Jd D x'd d k dG A (x-,x'h ^Af_(x-x'h k E (x') 


3 
4 


4b.II.II' 


- ! l^^ D ^^o^;x')7 fc ^A c + f _(x-x')7°2o(x / ) 


^8 


4b.II.II 


^^/ C / D x'<9 fc ^ix;x')7 fc ^A; f _(x-x')7^ (x') 


3 
8 


4b.II'.I 


- I ^Jd D x'd d k SG A (x ] x') 1 °^A'f_(x-x') 1 k E (x') 


+1 


4b.II.I 


-^^Id D x'd k d e SG A (x;x') 1 k ^Af_(x~x') 1 e E (x') 


+1 


4b.I.ir 


I ^Jd D x , d k d f dG A (x-x / ) 1 k ^iAf_(x-x / ) 1 E Q (x') 


1 

8 


4b.I.II 


- I( ^Id D x'd k d e SG A (x ] x'h k ^A c _l_(x-x') 1 e E (x') 


+ 1 


4b.I.I.A 


-^Jd D x'd k d k dG A (x-,x')Y^A c l_{x-x')Y^o{x') 


~32 


4b.I.I.BC 


-\i^Id D x'd k d e 6G A (x ] x') 1 k ^A^x-x') 1 ^o(x') 


1 

32 


Total 


-ifjd D x'd a d , ^dG A (x ] x')[ i ^%J j ^iAf_^ 


5 

16 



Table 6: The full contribution from each (46) term consists of its numerical 
coefficient times ll f & H 2 Haj E (x). 



Now extend the range in the first integral and take D = 4 in the second, 
fdtt%- 3 (l-t)- D+3 

a 

= f 1 dtt%- 2 (l-t)- D+3 + Cdt- + 0(-,D-4) = - + 0(1) . (145) 

a 

This gives, 

1 1 vi 3 3 

I 2 — > — e x 47r x aH x - x i—e x — = — 7r 2 x iaH . (146) 

2 e 2 8e 8 

Substituting our results for 7i_ 4 into expressions (I126H129P gives the en- 
tries for 2b.II'.I, 2b.II.I, 2b.I.II' and 2b.I.II in Table EJ Combining the totals 
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Term 


Coef. 


6 





7 


3 


8 






Table 7: The full contribution from each term consists of its numerical coef- 



ficient times -H"q7°H 



: {x). 



from Tables [3H7| gives a result in perfect agreement with our explicit compu- 
tation (HOD , 



k 2 n 



{i$* 2 (x),tf(k,s)} 

f 1 1 



k 2 H 2 
16tt 2 



16 16 



iHaj°E (x) 

k 2 H 2 
16n 2 



iHa~f°Z (x) x — . (147) 
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6 Discussion 

We have taken a major step in developing a technique to sum the series 
of leading infrared logarithms of inflationary quantum gravity. Our tech- 
nique was to employ a previous explicit computation [3l H] as "data" in the 
search for a simple operator formalism for reproducing the leading infrared 
logarithms. We found that only gravitons with the A-type polarization con- 
tribute, and only a single term fl85l) in the infinite series expansion of their 
propagator matters. We do not yet know if our new rule (I86H87P reproduces 
the leading logarithms of other quantities, or if it continues to work beyond 
one loop for the fermion effective mode function. 

One can easily see that the old rule ([318]) fails to reproduce the leading 
logarithms of quantum gravity. For example, consider the term l.b.II of 
Table El 



lb.II = 



IK 



2(D - 3) 
k 2 H 2 



d D x'd k iSA? + (x; x')?? G c f(x-x')j k E (x') , (148) 
1 



16tt 2 



iaHj Eq(x) x 



(149) 
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The old replacement (J7J) corresponds to substituting a purely temporal func- 



tion for i5A A Ax\ x 



iSAtix; x') — + H ° I r(jP n ^ 2 In 
+A ' (4tt)t r(f) 



mm (a, a 



(150) 



Of course the spatial derivative of this would give zero for lb. II! 

There are also problems with the closely related contribution from 3b. II, 

3b.II = ~ 2{ ^_ 3) / d D x'd k dG A (x; x')^tA?ix-x'h k ~ (x') , (151) 
The old replacement (jSJ) corresponds to the substitution, 

This gives a nonzero result, but not the right one, 

m 2 r(f + 1) „ . . n diie^ 

7°^o(x) / , , (154) 



2n^H D-3 ' ^ J J-^a' D ^A v D +^ 
l lJL l aH 1 Eo(x) x -4. (155) 

l07T Z 



In fact the old rule ([318]) does not give correct results for any of the thirty 
distinct nonlocal contributions of Tables [3]l6] The failure of this rule — which 
works for models without derivative couplings [TUJ [231 121] — deserves com- 
ment. Massless, minimally coupled scalars and gravitons are active fields. 
In order to produce infrared logarithms a theory must possess interactions 
involving at least one undifferentiated active field. However, there is a hier- 
archy of increasingly complicated ways in which this can happen: 

1. The theory may involve only undifferentiated active fields; 

2. The theory may involve active and passive fields with non-derivative 
interactions; and 
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3. The theory may involve differentiated active fields, with or without 
passive fields. 

The relation of our new rule (I86H87I) to the old rule ([HIED can be understood 
by considering how the expectation value of a given term in the free field 
expansion of some operator attains leading logarithm order in each case. 

When only undifferentiated active fields are present, each pair of free fields 
and each vertex integration must contribute to an infrared logarithm [24]. 
Therefore only the infrared part of the free field mode sum matters and one 
can effect this truncation at the level of the Yang-Feldman equations. This is 
the case solved by Starobinskii and Yokoyama [23J . At the level of expectation 
values of the free field expansion it corresponds to the replacements (J7HED 
with D = 4 because there are no ultraviolet divergences at leading logarithm 
order. 

When passive fields are present, but the active fields are not differen- 
tiated, reaching leading logarithm order still requires every active field or 
active Green's function to contribute to an infrared logarithm. Passive fields 
cannot be infrared truncated because the order [ln(a)]° contributions they 
make derive from all parts of the free field mode sum. However, precisely 
because passive fields cannot produce infrared logarithms we can perform pas- 
sive vertex integrations without accounting for the spacetime dependence of 
active fields or Green's functions. This amounts to integrating out the passive 
fields and then evaluating the resulting, nonlocal effective action assuming 
the active fields are constant — which defines the effective potential. At the 
level of expectation values of the free field expansion it corresponds to the 
replacements ([HIE]) with dimensional regularization on because there can be 
ultraviolet divergences at leading logarithm order. 

The situation is vastly more complicated when differentiated active fields 
are present. In this case the vertex integration of a differentiated active 
field propagator or Green's function can produce an infrared logarithm, even 
though the integrand contains no logarithm. For the fermion wave func- 
tion, every infrared logarithm arises in this fashion. We cannot ignore differ- 
entiated active fields because they can still contribute infrared logarithms. 
Neither can we ignore their spacetime dependence in performing vertex in- 
tegrations, and we must retain dimensional regularization in order to define 
these integrals. In view of this it seems doubtful that any infrared truncated 
formalism can correctly represent the theory, even at leading logarithm order. 
The replacements (I86H87P of our new rule seem to represent the appropriate 
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generalizations of the old rule (JTHHD to this more singular environment. 

In addition to showing that the new rule works, our analysis provides a 
deeper understanding of why the fermion mode function acquires a secular 
enhancement whereas the scalar mode function does not [36J. The reason is 
spin. At late times the kinetic energies of all quanta redshift to zero. This is 
why we could neglect the d' j a E (x') contributions from terms (la), (2a), (3a) 
and (4a) of Table [2J A massless scalar interacts with gravity only through 
its kinetic energy. Inflationary particle production immerses such a scalar 
in a sea of infrared gravitons but they do little because the interaction is so 
weak. In contrast, a massless fermion possesses an additional gravitational 
interaction through its spin, which does not redshift. That is why we found 
leading order contributions from the 7 P J tTj aS (x / ) terms of (16), (26), (36) and 
(46) on Table H 

Gravitons also have spin and it is natural to wonder what the sea of 
infrared gravtions does to itself. One could answer this by using the known 
one loop graviton self-energy [11] to correct the graviton mode functions, 
just as we have done for fermions. It would also be interesting to understand 
in this way the null result that has been obtained at one loop order for the 
graviton 1-point function |37J. In particular, can the spin-spin interaction 
lead to significant quantum gravitational back-reaction? 
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